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ABSTRACT 


We investigate whether the sky rate of Fast Radio Bursts depends on Galactic latitude using the first 
catalog of Fast Radio Bursts (FRBs) detected by the Canadian Hydrogen Intensity Mapping Experi- 
ment Fast Radio Burst (CHIME/FRB) Project. We first select CHIME/FRB events above a specified 
sensitivity threshold in consideration of the radiometer equation, and then compare these detections 
with the expected cumulative time-weighted exposure using Anderson-Darling and Kolmogrov-Smirnov 
tests. These tests are consistent with the null hypothesis that FRBs are distributed without Galactic 
latitude dependence (p-values distributed from 0.05 to 0.99, depending on completeness threshold). Ad- 
ditionally, we compare rates in intermediate latitudes (|b| < 15°) with high latitudes using a Bayesian 
framework, treating the question as a biased coin-flipping experiment- again for a range of completeness 
thresholds. In these tests the isotropic model is significantly favored (Bayes factors ranging from 3.3 to 
14.2). Our results are consistent with FRBs originating from an isotropic population of extragalactic 


Sources. 


Keywords: Radio transient sources (2008); High energy astrophysics (739) 


1. INTRODUCTION 


Fast Radio Bursts (FRBs) are us—ms duration radio 
bursts with dispersion measures (DMs) consistent with 
origins far outside our Milky Way galaxy. Thirteen FRB 
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sky localizations (Tendulkar et al. 2017; Chatterjee et al. 
2017; Bannister et al. 2019; Ravi et al. 2019; Macquart 
et al. 2020; Marcote et al. 2020; Heintz et al. 2020) have 
identified host galaxies, most at cosmological distances, 
but which are of diverse types that have not yet clari- 
fied the nature of these transient sources. Some FRBs 
repeat (Spitler et al. 2016; CHIME/FRB Collaboration 
et al. 2019a,b; Fonseca et al. 2020; Kumar et al. 2019), 
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and repeating sources exhibit bursts having properties 
that differ temporally and spectrally from those of ap- 
parent non-repeaters (Spitler et al. 2016; CHIME/FRB 
Collaboration et al. 2019b; Hessels et al. 2019; Fonseca 
et al. 2020; Pleunis et. al. 2021). Hence, there may be 
more than one class of bursting radio sources comprising 
FRBs. 

One puzzling observational feature of FRBs early on 
in their study was evidence for a Galactic latitude de- 
pendence of their sky distribution. Petroff et al. (2014) 
reported that FRBs avoid the Galactic plane at a 99% 
confidence-level based on four events detected at high 
latitudes and a lack of detections in an intermediate 
latitude FRB survey (—15° < b < 15°), with both 
surveys done using the 64-m Murriyang telescope! at 
the Parkes Observatory. Further observations report- 
ing an additional two bursts at high latitude supported 
this claim, arguing for a ~ 3c dearth at low latitudes 
(Burke-Spolaor & Bannister 2014). Champion et al. 
(2016) found, following the detection of an additional 
five events, that the intermediate and high latitude rates 
are inconsistent with 97.596 confidence. However, Con- 
nor et al. (2016) used the same data set but different 
statistical techniques to argue against any latitude de- 
pendence. Rane et al. (2016) also argued against a de- 
pendence, attributing the perceived deficit to initially 
overestimated FRB rates. Most recently, based on a 
sample of fifteen FRBs, Bhandari et al. (2018) concluded 
that any Galactic latitude effect has low (< 2c) signif- 
icance, and that past claims could be explained by ei- 
ther small samples of FRBs, hence poor statistics, or 
inaccuracies in the modeling of Galactic dispersion and 
scattering effects. 

Meanwhile, Macquart & Johnston (2015) suggested 
that diffractive scintillation could effectively boost the 
apparent number of high-latitude sources even for an un- 
derlying isotropic distribution. However this requires a 
very steep FRB flux density distribution. The scarcity of 
Arecibo FRB discoveries provided evidence against this 
distribution, given the Arecibo telescope’s high sensitiv- 
ity (Scholz et al. 2016). While unlikely given the high 
DMs of the Parkes detections, another generic reason 
for a low-level Galactic latitude dependence could be 
the presence of a relatively small fraction of Galactic 
sources within the FRB population at high latitudes, 
on top of a largely isotropic distribution, with the for- 
mer misidentified as extragalactic due to inaccuracies in 
Galactic DM-distribution models. 


1 Formerly known as the Parkes Radio Telescope. 


The Canadian Hydrogen Intensity Mapping Exper- 
iment (CHIME?) Fast Radio Burst project (hereafter 
CHIME/FRB; CHIME/FRB Collaboration et al. 2018) 
is well suited to searching for Galactic latitude de- 
pendence of the FRB sky distribution because of the 
system's high event rate and uniform, repeated cover- 
age of the sky northward of declination —11?. Oper- 
ating since mid-2018 in the 400-800 MHz band, the 
first CHIME/FRB catalog (CHIME/FRB Collaboration 
2021), hereafter Catalog 1, reports on 535 events de- 
tected with a single search pipeline, hence uniform se- 
lection biases. Here we report on our search for any 
Galactic latitude dependence of the FRB sky distribu- 
tion using the Catalog 1 sample. 

The paper's outline is as follows. In 82 the observa- 
tions are discussed, namely the exposure maps and the 
sample of CHIME/FRB events included in the analy- 
sis. Sensitivity correction methods are presented in 83. 
Results from two classes of statistical tests are given in 
$4 and outcomes are discussed in 85. Conclusions are 
provided in 86. 


2. OBSERVATIONS 


'The primary observational inputs taken from Catalog 
1 are the exposure map, source positions, and detec- 
tion signal-to-noise ratios (S/N; used for completeness 
cuts). Exposure maps are provided in a HEALPix? for- 
mat (Górski et al. 2005) and span from 2018 August 
28 to 2019 July 1. Notably, the integration range does 
not cover all bursts in Catalog 1, as the reported ex- 
posure is meant to capture relatively uniform observing 
conditions (e.g. singular beam configuration and sta- 
ble sensitivity). The standard of uniformity is crucial 
for this analysis, so we excise events accordingly. The 
first thirteen bursts (CHIME/FRB Collaboration et al. 
2019c) were detected during a precommissioning phase 
with a considerably different beam configuration. A fur- 
ther 26 events occurred on days with either abnormally 
low sensitivity, days with disruptive hardware/software 
upgrades, or days with significant pipeline issues. Three 
events are removed for being detected in the far side- 
lobes of the telescope, as obtaining their positions re- 
quires methods beyond the scope of Catalog 1. Finally, 
we include only a single burst from each of the eighteen 
repeating FRB sources in Catalog 1, taking the event 
with the highest S/N (this choice is to accommodate 
subsequent completeness cuts, see 83.5). The final sam- 
ple used in this paper includes 453 bursts. The exposure 


? https:/ /chime-experiment.ca 


3 https:/ /healpix.sourceforge.io 
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products provided by Catalog 1 have been defined at the 
FWHM of the synthesized beams at 600 MHz, which ne- 
cessitates a high map resolution to capture gaps in the 
exposure between beams. Such a resolution is not re- 
quired for this work, so we downsample the maps to 
have pixel area ~ 47 square arcminutes and use the 
same resolution when generating the sensitivity maps 
used for completeness cuts (see §3). The beam con- 
figuration for this exposure spans 120° North/South in 
elevation angle centered at zenith, so sources with a dec- 
lination greater than ~ + 70° are observed twice a day. 
Because of this, separate exposure maps exist for the up- 
per and lower transits. Cumulative exposure and source 
positions along Galactic latitude and declination coor- 
dinates are shown in Figure 1 (with only the cuts de- 
scribed in this section applied to the data). The sum 
total of the exposure is ~6.6 x 10° hours: deg?. Note 
that the analysis presented in this paper uses the S/N 
as the relevant quantity rather than fluence when mak- 
ing completeness cuts. The current localization uncer- 
tainties lead to large uncertainties in the fluence mea- 
surements, which are treated as lower limits as bursts 
are assumed to occur in the center of the formed beams 
for fluence calculations. As part of Catalog 1, a fluence 
completeness threshold is estimated for each burst via 
Monte Carlo simulation. These methods aim to cap- 
ture a range of sensitivity variation and are described in 
detail by Josephy et al. (2019). The reported complete- 
ness thresholds are 95% one-sided confidence intervals, 
and the median value of these thresholds is ~5 Jy - ms. 
The triggering threshold of the real time pipeline is a 
S/N of 9, and we adopt this value in the subsequent 
analysis for burst detectability. As for the DM trigger- 
ing criteria, CHIME/FRB will label an event as extra- 
galactic when the DM is larger than either of the two 
main Galactic free-electron models estimates (NE2001: 
Cordes & Lazio 2002; YMW16: Yao et al. 2017), using 
the formed beam center as the line-of-sight. Note that 
this is different from initial plans (CHIME/FRB Collab- 
oration et al. 2018). The change was made early in the 
precommissioning phase and has the advantage of being 
simpler, more generous, and without reference to mixed 
uncertainty estimations. 


3. SENSITIVITY MAPS 


From Figure 1, a naive comparison of the cumulative 
CHIME/FRB events against exposure does not quali- 
tatively suggest a Galactic latitude dependence. How- 
ever, there is a noticeable deviation from isotropy when 
declination is taken as the coordinate of interest. This 
deviation is not a surprise given the significant sensi- 
tivity fall-off in the primary beam with zenith angle 
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Figure 1. Cumulative exposure (dashed) and cumulative 
source positions (solid) as a function of absolute Galactic 
latitude (black) and as a function of declination (blue). Note 
the slope increase above declination ~70° is due to doubled 
exposure from the secondary "lower" transit. Only the cuts 
described in 82 have been applied to the original Catalog 1 
sample of bursts. 


(see Figure 2). While there is an interesting astrophys- 
ical motivation for examining Galactic latitude rate de- 
pendence (e.g., evaluating Galactic free-electron density 
models), a rate dependence with declination must be in- 
strumental in the case of CHIME/FRB. Conducting the 
statistical tests against both Galactic latitude and dec- 
lination coordinates in parallel serves as an important 
cross-check that instrumental effects are not corrupting 
the main focus of this study. To properly compare ex- 
posure and detections, we aim to filter out events be- 
low a completeness threshold. In the context of this 
work, we say a burst is “above completeness” if we ex- 
pect it to be reliably detected throughout the given field 
of view. We can make this comparison with a relative 
sensitivity map, which we model using a reduced version 
of the single-pulse form of the radiometer equation (e.g., 
Cordes & McLaughlin 2003): 


Tsys vU, (1) 
Gv Av Wi 


Here Smin is the minimum detectable flux density, 
Tsys is the system temperature, G is the telescope gain, 
Av is the observing bandwidth, w, is the broadened 
pulse width, and w; is the intrinsic pulse width. Note 
that factors common to the CHIME/FRB instrument 
have been dropped here, namely the number of polar- 
izations, S/N threshold, and scalar efficiency (i.e., digiti- 
zation loss). As the scope of this work is limited to inves- 
tigating a differential rate in the CHIME/FRB sample, 
we can afford to operate in relative terms, which simpli- 
fies the analysis substantially. The system temperature, 


Smin e 
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forward gain, effective bandwidth, and pulse-broadening 
terms are treated separately. Each term leads to a rela- 
tive sensitivity map in HEALPix format and a compos- 
ite map is formed through simple multiplicative combi- 
nation. Just as separate exposure maps exist for the up- 
per and lower transit, we produce two sensitivity maps 
to capture the double-valued gain. The components of 
the composite map are described below. 


3.1. System Temperature 


Sensitivity is linearly proportional to the system tem- 
perature, which we take to be Tsys = Tsky + 50K using 
the receiver temperature from CHIME/FRB Collabora- 
tion et al. (2018). We estimate sky temperature starting 
with the destriped, but not desourced, Haslam all-sky 
continuum map at 408 MHz (Remazeilles et al. 2015). 
The map is dominated by synchrotron radiation, which 
is well described by a power-law frequency dependence. 
We use the spectral index map from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP; Bennett et al. 2003) 
in conjunction with the Haslam map to compute a band- 
averaged sky temperature. The mean temperatures for 
intermediate (|b| < 15?) and high Galactic latitudes are 
21 K and 10 K respectively, resulting in a ~18% mean 
difference in sensitivity. Note that, for the frequencies 
(1.22—1.52 GHz) and field of view relevant to Petroff 
et al. (2014), the differential sensitivity due to sky tem- 
perature was much less significant. 


3.2. Gain 


'The sensitivity varies dramatically across the field of 
view. The most apparent effect is a roll-off with increas- 
ing zenith angle. We use a composite beam model that 
includes a frequency-dependent description of both the 
primary beam and the synthesized beam pattern. The 
synthesized component is described by Ng et al. (2017), 
while the most recent description of the primary com- 
ponent is given by CHIME/FRB Collaboration (2021). 
'The model computes a normalized sensitivity after spec- 
ifying a frequency and position. To generate a senstivity 
map, we compute a beam and band averaged sensitivity 
for a fine gridding along the meridian (see Figure 2). Av- 
eraging for the formed beams is done within the FWHM 
at 600 MHz. 


3.3. Bandwidth 


Minimum detectable fluence is inversely proportional 
to the square root of the bandwidth. The total avail- 
able band fluctuates over time due to correlator clus- 
ter health (Denman et al. 2020), and is further re- 
duced by realtime masking of radio frequency interfer- 
ence (RFI). To quantify any spatial bandwidth depen- 
dence, we recorded metrics of the masking fraction as a 


Declination (deg.) 
-10 0 10 20 30 40 50 60 70 80 90 80 70 


-60  -40 -20 0 20 40 60 
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Figure 2. Normalized band-averaged gain (intensity beam 
response) along meridian. Sensitivity obtained with a com- 
posite beam model which includes primary and formed beam 
effects. Each point represents an average sensitivity within 
a synthesized beam's 600- MHz FWHM. The dashed portion 
of the profile shows the lower transit. The sensitivity spike 
at zenith is a consequence of the beamforming technique (Ng 
et al. 2017). 


function of frequency (downsampled to 1024 channels), 
on a per-beam basis, downsampled from 1 ms to 60 s. 
Since we are building towards a singular sensitivity map 
that reflects bandwidth losses due to persistent RFI, we 
average across the four East- West beam columns and av- 
erage in time across several metric collection runs from 
2020 March 22 to 2020 July 30 (see Figure 3). Correlator 
related fluctuations are included in these metrics, but 
affect all formed beams equally and are subdominant 
to RFI masking, so we do not pursue an independent 
characterization of the correlator health over time. The 
resulting mask has a complicated spatial structure with 
an increase in masking towards the horizon. The differ- 
ential effective bandwidth amounts to a relative sensi- 
tivity loss of roughly 10%. The ability to stream and 
aggregate the RFI metrics was implemented after the 
2019 July 1 cutoff date for the exposure. To investigate 
whether or not the RFI environment had changed signif- 
icantly, we visually compared the time-averaged masks 
with the per-event masks, which are saved along with 
intensity data for each event. The enhanced masking 
from 450 MHz to 600 MHz near the horizons and the 
midband hot-spot that is seen in Fig. 3 are also appar- 
ent in the per-event masks after averaging event masks 
by beam row, suggesting that direction dependent RFI 
conditions have not changed considerably. 


3.4. Pulse Broadening 


Unlike the preceding sensitivity corrections, effects 
due to pulse-broadening are difficult to apply uniformly 
as they vary on a pulse-to-pulse basis. A narrow pulse 
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Figure 3. RFI masking fraction as a function of zenith angle 
and frequency. Masking metrics are averaged over the four 
beams in the East-West direction. The top panel shows the 
square-root of the normalized bandwidth, which translates 
to a relative sensitivity. 


originating relatively nearby will be scatter-broadened 
by the Galaxy to a larger fractional degree as compared 
to a distant and intrinsically broad pulse. We take the 
broadened width as the quadrature sum 


uy = yw? + tamp + Bhan + 72 (2) 
where w; is the intrinsic width, tsamp is the sampling 
time (0.983 ms), Tsc is the scattering timescale, and tehan 


is the intra-channel dispersion smearing, 


AVehan po DM 
tehan = 8. —— — — —M is 
n iras ( MHz ) (aE) (= . ~) (3) 


Here Avenan is the channel bandwidth for CHIME/FRB 
(24.4 kHz) and v is the channel frequency. For refer- 
ence, the intra-channel smearing at 600 MHz for DM 
= 500 pc cm? is ~0.5 ms for CHIME/FRB. Galactic 
scattering timescales are estimated with electron den- 
sity models (NE2001: Cordes & Lazio 2002; YMW16: 
Yao et al. 2017) and we use a Tse x v^? relationship 
to appropriately scale from 1 GHz (Bhat et al. 2004). 
We simplify Galactic scattering as originating from a 
screen located 25 kpc away and include the geometrical 
lever-arm effect (Lorimer et al. 2013) 


Tsc — Tmas : 4f (1 — f); f= zh (4) 
L 

where Tmar is the maximum scattering timescale 
achieved while integrating the models out to 25 kpc 
and Dry, is the luminosity distance to the FRB source. 
Altogether, a relative sensitivity map due to pulse- 
broadening depends on w;, DM, Dr, and v. Since such 
a map is burst dependent, we construct an averaged 


map using Monte Carlo methods, taking a log-normal 
fit of the intrinsic widths from Catalog 1 to specify a 
width distribution. For a distance we take 100 Mpc for 
all bursts, despite most FRBs being located at much 
further distances. ‘This simplified treatment is moti- 
vated by the asymptotic behaviour of the lever-arm ef- 
fect. Even with a conservative placement of all FRBs 
at 100 Mpc, the final averaged sensitivity map shows 
that pulse-broadening effects are subdominant to other 
sensitivity factors relevant to CHIME/FRB. The mean 
relative sensitivity within intermediate latitudes is 9996 
and effectively 10096 at high latitudes. Consequently, we 
should not over-interpret a lack of latitude dependence 
as an outright validation of the Galactic free-electron 
models. If we consider the other effects to be a fore- 
ground to pulse-broadening effects, we see how impor- 
tant modeling the beam, sky temperature, and band- 
width are. To this end, it will be worthwhile repeating 
the analysis for future CHIME/FRB catalogs with up- 
dated components (e.g., a sky temperature map mea- 
sured directly by the CHIME cosmology experiment). 
Note that we have not considered any host-related scat- 
tering here, as those should not show latitude depen- 
dence and can be absorbed into the intrinsic width term. 


3.5. Completeness Cuts 


With HEALPix maps representing relative sensitivi- 
ties due to differential sky temperature, effective band- 
width, forward gain, and pulse-broadening, a composite 
sensitivity map can be made with a simple multiplica- 
tive combination. This map is normalized to its most 
sensitive position (see Figure 4). Note that, beyond ex- 
cising bursts from days with anomalous sensitivity (see 
$2), we do not incorporate any day-to-day sensitivity 
corrections when filtering events for completeness. We 
say a burst is above completeness when we expect it to 
be detectable throughout the field of view, so the pro- 
cess of forming a completeness cut requires the sky area 
to be defined apriori. Rather than specifying the field of 
view explicitly, we make our completeness cuts by first 
designating a minimum sensitivity threshold (the cut). 
Locations in the map with a sensitivity lower than this 
cut are excluded from the field of view (exposure pix- 
els are set to zero). An FRB with non-zero exposure 
is then considered above completeness if the following 
inequality is true 


x Sensitivity cut 
Sensitivity at FRB 


S/Nprg > 9. (5) 
A single choice for the sensitivity cut is somewhat arbi- 
trary, although natural choices might be at half of peak 
sensitivity, whichever cut maximizes the sample of FRBs 
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above completeness, or some value that balances expo- 
sure with sample size (see Figure 5). There are two com- 
peting effects that dictate the final sample size. First, 
as the sensitivity cut is increased, the field of view is 
decreased and the number of FRBs with non-zero ex- 
posure shrinks. Secondly, as the cut increases, the S/N 
scaling becomes less extreme and a larger fraction of the 
FRBs with non-zero exposure remain detectable. As the 
choice of cut is arbitrary, we repeat the statistical tests 
for latitude dependence for a range of sensitivity cuts. 


Galactic coordinates; Mollweide projection 


Latitude 


-75? 


Longitude 


LS 
00 02 04 06 08 1.0 
Normalized Sensitivity 


Figure 4. Normalized sensitivity for the CHIME field 
of view, including contributions from sky temperature, in- 
tensity beam response, masking fraction, DM intra-channel 
smearing, and scattering. Normalization is done with respect 
to the peak sensitivity. The map has been rotated 123° in 
Galactic longitude to horizontally center the North Celestial 
Pole. NE2001 is used here for pulse-broadening terms, with 
YMW/106 giving a qualitatively similar map. White dots show 
the location of the FRBs used in this analysis. Dashed black 
line shows the extent of the lower transit, where exposure 
and sensitivity is double-valued (lower transit sensitivity not 
shown). The concentric circular patterns are due to beam 
sensitivity (Fig. 2) and RFI masking (Fig. 3). 


4. STATISTICAL TESTS 
4.1. Cumulative Distribution Tests 


Assuming sensitivity corrections have been properly 
handled, the detections above completeness should track 
exposure if the FRB rate is isotropic across the sky. Sta- 
tistical tests that ask whether an empirical distribution 
function (EDF) is consistent with a cumulative distribu- 
tion function (CDF) are natural to apply here. In par- 
ticular, we use non-parametric one-sample Kolmogorov- 
Smirnov (K-S) and Anderson-Darling (A-D) tests. A-D 
tests are often considered preferable for certain distri- 
butions, where deviations in the tails are not realized 
by the difference statistic of a K-S test. Since Bhandari 
et al. (2018) concluded that a departure from isotropy 
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Figure 5. Number of FRBs above completeness (black 
solid) and sky area (orange dashed) as a function of sen- 
sitivity cutoff. As the threshold is raised, the sky area and 
the number of candidate FRBs for inclusion is reduced. Like- 
wise, fewer candidates are filtered out for having insufficient 
S/N to survive a translation to a location with the threshold 
sensitivity. Whether NE2001 or YMWI16 is used to calculate 
pulse-broadening effects, the samples sizes are nearly identi- 
cal across all sensitivity cuts. 


was not significant via a K-S test, we include it here 
for comparison. The exposure CDF and detection EDF 
that both tests rely on are obtained as follows. After 
choosing a coordinate along which the comparison is to 
be made (e.g., declination or Galactic latitude), the ref- 
erence CDF is generated from the exposure maps by 
weighting the coordinate value at each pixel by its ex- 
posure value. Pixels with a relative sensitivity below 
the chosen threshold are assigned a weight of zero. The 
EDF follows from the positions of the events above com- 
pleteness. Examples of sensitivity corrected cumulative 
exposure and detection distributions are shown in Fig- 
ure 6 for both absolute Galactic latitude and declina- 
tion, while test results for a range of sensitivity cuts are 
shown in Figure 7. Given the multitude of tests, their 
interpretation is not so straight-forward. Trial factors 
or look-elsewhere effects are difficult to quantify as tests 
performed using different sensitivity cuts are not inde- 
pendent (for similar thresholds, the exposure CDF and 
detection EDF are also similar, hence similar p-values). 
A relevant statistic is the harmonic mean p-value (HMP; 
Wilson 2019), which controls for familywise error rate 
and, importantly, is robust to inter-dependent samples 
of p-values. For N equally weighted tests, the HMP is 


N 
3 5 1/pi 


The HMP should only be interpreted as a false-alarm 
rate in the usual way for small N and small p. Adjusted 


(6) 


p= 


GALACTIC LATITUDE DEPENDENCE OF FRB Sky DISTRIBUTION 7 


critical values are obtained with with the CRAN pack- 
age harmonicmeanp*. Considering the 80 tests shown 
in Figure 7, the adjusted critical value for a 0.05 signif- 
icance threshold is 0.037. The HMP for the A-D tests 
is 0.56, so we conclude the ensemble of tests does not 
support a Galactic latitude dependence. 


4.2. Bayesian Coin Flipping Tests 


Connor et al. (2016) equate the problem of determin- 
ing whether a differential rate exists to a biased coin 
flipping experiment and provide a Bayesian framework 
to test for isotropy, which we follow here. The sky is 
first partitioned in two, separating high-latitude and 
low-latitude regions about |b| = 15°. For an isotropic 
sky distribution (Model 1), the probability of detect- 
ing a random burst in either region is dictated by their 
relative exposure. With a as the ratio of high-latitude 
exposure to low-latitude exposure, the probability of a 
random burst appearing in the high-latitude region is 
p = a/(a+1). The probability of detecting M out 
of a total K events in the upper region then follows a 
binomial distribution 


PMK p) -(Sa-»". m 


'To conduct a Bayesian model selection, Connor et al. 
(2016) specify an alternative possibility (Model 2), 
where p is not defined by the relative exposure between 
regions, but is instead an unknown free parameter with 
a flat prior over 0 — 1. Assigning equal weights to each 
model, and marginalizing over the possible values of p in 
the alternative model, the ratio of posterior probabilities 
(Bayes factor) for the two models is 


P(M|Model 1, kK) - K\ m "T 
P(M|Model 2, K) - T (1—p1)* ^. (8) 


Once again we perform this model selection test for 
a range of sensitivity cuts, where each cut defines the 
values of K, M, and a (which gives pı). Resulting Bayes 
factors along with the exposure ratio a and detection 
ratio K/(M — K) are shown in Figure 8. To interpret 
the significance of the resulting Bayes factors, we use 
the scale suggested by Kass & Raftery (1995), where 
the evidence for Model 1 over Model 2 is “strong” when 
the Bayes factor is greater than 10, “substantial” when 
between 10!/? and 10, and “barely worth mentioning” 
between 1 and 10!/2. We find the evidence for isotropy 
is strong for roughly half of the sensitivity cuts, and 
substantial for the remaining half. 


4 https://CRAN.R-project.org/package=harmonicmeanp 


4.3. Cross-checks 


Given the large size of the initial sample of 
CHIME/FRB events, we are afforded the luxury of per- 
forming a number of additional cross-checks that would 
be difficult with previous FRB samples. We can check 
the robustness or stability of results by jackknife resam- 
pling, and we can partition the sample with respect 
to a variety of parameters to check for possible lati- 
tude deficits within subsections of the parameter-space. 
For example, intrinsically narrow bursts may suffer from 
Galactic scatter-broadening or higher frequency bursts 
may be suppressed through scintillation effects. Test re- 
sults for sample splits are summarized in Table 2. Of 
note here is the upper excess DM split (DMx > me- 
dian DM x), where the HMP is below the adjusted crit- 
ical value for 0.05 significance. When compared to ini- 
tial claims, the significance is somewhat marginal (false- 
alarm rate of ~2%). It is not clear why events with 
higher excess DM near the plane are under-represented 
in our sample and it will be interesting to revisit this 
analysis with future CHIME/FRB catalogs to see if the 
effect persists. 

We have also investigated splitting the sample into 
one-offs and repeaters. The results are effectively the 
same with respect to the A-D tests, while for the Bayes 
tests, the previously strong evidence for isotropy is di- 
minished to being “barely worth mentioning". How- 
ever, we caution that the sample size of repeaters used 
in this analysis is quite small in comparison to the one- 
offs. With only order 10 sources surviving the sensitiv- 
ity cuts, the current statistical power is quite border- 
line. CHIME/FRB will continue to expand the pool of 
repeaters, and this analysis should be revisited as the 
population grows. 

A-D test results along Galactic and celestial coordi- 
nates are summarized in Table 1, including tests against 
both absolute and signed Galactic latitude. Elsewhere, 
we have chosen to focus on absolute latitude by default, 
under the assumption that there is no sign dependence, 
and that a more dense sampling of the parameter space 
is statistically preferable. In both cases the HMP is 
unremarkable, suggesting our symmetry assumption is 
valid. Deviations from isotropy with respect to declina- 
tion are indicative of unmodeled instrumental effects, so 
the excess of events just below declination 6 ~ 30° is 
cause for concern (see Figure 6, bottom row). For decli- 
nation, the A-D p-value is minimal (0.01) at a sensitiv- 
ity cut of 0.32. The HMP is above the adjusted critical 
value for a 0.05 significance, so we say the ensemble of 
declination tests does not reject isotropy. Nevertheless, 
we investigate the worst-case sensitivity cut by taking 
half of the post-cut events randomly (without replace- 
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Figure 6. Distributions of detection latitudes (empirical, solid lines) and exposure (cumulative, dashed lines). Absolute Galactic 
latitude is shown on the top row and declination on the bottom. Each column corresponds to a different sensitivity cut and 
filtered sample of detections. Shaded regions are a 2e band derived from an ensemble of simulated detections, where each mock 
sample includes the same number of bursts as the observed sample and is generated from the exposure CDF. 
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Figure 7. Anderson-Darling (solid lines) and Kol- 


mogorov-Smirnov (dashed lines) test results for a range of 
sensitivity cuts. Tests compare cumulative exposure to cu- 
mulative detections (above completeness) as a function of 
absolute Galactic latitude. 


ment) and repeating the A-D test 1000 times. The mean 
p-value of these jackknife resampling tests is 0.12, sug- 
gesting the deviation from isotropy is marginal. Further- 
more, splitting the sample on arrival time, the latter half 
does not show a statistically significant excess. It is not 
clear what would explain the excess- events in that decli- 
nation range are not clustered in right ascension or DM, 
so a missed repeater does not work as an explanation for 
the excess. An increased sensitivity may be expected for 


declinations close to the sources responsible for the daily 
system-wide gain calibration. 'These sources: Cygnus 
A (ó ~ +41°), Cassiopeia A (ô ~ +59°), Taurus A 
(6 ~ +22°), and Virgo A (6 ~ +12°), are cycled depend- 
ing on when their transits occur. During the first half of 
Catalog 1 Taurus A was not used, so an increased sen- 
sitivity in that declination regime is not expected. Re- 
peating the A-D tests for Galactic latitude using only 
the latter half, the p-values remain insignificant. As 
these methods are sensitive to both scarcity and ex- 
cess, it is reasonable to investigate questions of possible 
Galactic contamination. A concrete limit would be rele- 
vant to evaluations of NE2001 and YMW16, and would 
be a useful figure of merit for any low DM population 
inferences of Catalog 1. Fully quantifying what frac- 
tion of the sample could be of Galactic origin is beyond 
the scope of this paper, as to do so accurately would 
require integrating the CHIME/FRB selection function 
with population synthesis models. However, working 
towards a rough approximation, we have considered a 
simple limiting case where selection effects are ignored, 
and the contaminating population is composed entirely 
of young isolated neutron stars. This scenario concen- 
trates potential sources to low Galactic latitude, which 
should have the strongest statistical impact. We first 
simulate a large population of fake sources following 
the prescription from Faucher-Giguére & Kaspi (2006). 
After projecting onto the sky, this population can be 
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Figure 8. Diagnostic plots for Bayesian biased coin flip tests 
performed for a variety of sensitivity cuts. Upper panel: the 
sky is partitioned into high and intermediate latitudes at 
|b] = 15° — ratios between these regions are shown for both 
the detections (solid line) and the exposure (dashed line). 
The shaded region is a 2e band derived from an ensemble of 
simulated “coin flips”, where each mock includes the same 
number of events as the observed sample and is generated 
from a binomial distribution with bias specified by the ex- 
posure ratio. Lower panel: Associated Bayes factor for each 
sensitivity cut (solid line), along with common significance 
thresholds (dashed lines), with “strong” evidence for isotropy 
at a Bayes factor of 10, “significant” at 10'/?, and “barely 
worth mentioning” at 1. 


randomly sampled using the exposure maps as weights. 
We then repeated the A-D tests shown in Figure 7, 
progressively increasing the number of fake sources we 
add to the real sample. We find the HMP passes the 
critical points equivalent to 5%, 1%, and 0.1% false- 
alarm rates, when the contaminated fraction of the final 
sample reaches roughly 7%, 8%, and 10% respectively. 
Again, this should only be taken as a crude estimate of 
how sensitive these methods are to varying degrees of 
Galactic contamination. A full analysis should incorpo- 
rate priors on allowable DM excess and impose selection 
effects on the synthesized population, which should in- 
clude evolved pulsars and RRATs. 


5. DISCUSSION 


Using 453 events from Catalog 1, we have searched 
for evidence of dependence of the FRB sky distribution 
on Galactic latitude, but have found none, as expected 
if FRBs are a cosmological population and selection ef- 
fects have been sufficiently accounted for. This is in con- 


A-D p-value 


min HMP max 


Galactic Longitude 0.40 0.64 0.97 
Latitude 0.02 0.28 0.77 
|Latitude| 0.05 0.18 0.99 

Right Ascension 0.11 048 0.98 
Declination 0.01 0.07 0.48 


Table 1. Min, harmonic mean (Eqn. 6), and max p-values 
from Anderson-Darling tests performed over a range of sen- 
sitivity cuts from 0.2 to 0.8. Each batch of tests compares 
cumulative detections with the cumulative exposure as a 
function of different coordinates. In none of the coordi- 
nate frames is the null hypothesis of isotropy consistently or 
strongly rejected. Note that the critical p-value of 0.05 for 
significance is adjusted to 0.037 for the HMP test. K-S test 
results are similar in all cases and are therefore not shown 
(see Figure 7). 


A-D p-value Bayes factor 


min HMP max min mean max 


No split 0.05 0.18 0.99 3.3 9.3 14.2 


One-offs 0.05 0.17 0.98 2.5 8.9 13.9 
Repeaters 0.09 0.25 0.89 0.4 2.3 3.5 


0.002 0.02 033 05 22 7.2 

DMx 
| 010 033 087 15 64 9.5 
„_ t O17 048 096 35 68 9.2 
“| 001 008 090 20 71 98 
œ 1 001 0.07 073 14 53 90 
“| 004 020 099 14 65 96 
, t 034 063 099 34 73 94 
D 0.0 0.07 0.98 23 56 10.7 


+ 


Table 2. Test results for (absolute) Galactic latitude depen- 
dence for different splits of the initial sample, where up/down 
arrows denote upper/lower splits, partitioned at the median 
value. The min/(harmonic)mean/max values summarize a 
range of tests with varying sensitivity threshold (from 0.2 to 
0.8 of peak sensitivity). DM, Tsc, and w; are relevant for 
pulse-broadening sensitivity corrections, while ve (peak fre- 
quency of spectral fit) has relevance for frequency dependent 
effects such as scintillation. Note that the critical p-value of 
0.05 for significance is adjusted to 0.037 for the HMP test. 
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trast to past reports (Petroff et al. 2014; Burke-Spolaor 
& Bannister 2014; Champion et al. 2016) which found 
avoidance of the Galactic Plane in surveys conducted at 
1.4 GHz, but which were based on far fewer events than 
considered here. 

Macquart & Johnston (2015) considered the previ- 
ously reported Plane avoidance and suggested that the 
effect was due to diffractive scintillation within the 
Milky Way, with the underlying FRB sky distribution 
being isotropic. The idea was that an effect associ- 
ated with Eddington bias artificially enhanced the high 
Galactic latitude rate, with the enhancement expected 
to be greater for a steeper FRB flux distribution. 

The FRB sky isotropy we report for Catalog 1 events 
may not, however, shed light on the effect considered by 
Macquart & Johnston (2015), who argued an enhance- 
ment at high latitudes was plausible if the scattering 
imposed by the Galaxy was substantially overestimated 
by electron density models, and extragalactic scattering 
was not so extreme (e.g., Tse < 1 ms at 1.4 GHz) as 
to alter the coherence properties and quench interstel- 
lar scintillation. With the basic requirement of having a 
decorrelation bandwidth Avg comparable or larger than 
roughly half the observing bandwidth, their argument 
was supported by the existence of several high-latitude 
pulsars with Avg ~ 300 MHz at 1.4 GHz in conjunction 
with the ~300 MHz bandwidth used in early FRB de- 
tections with the Murriyang telescope. While these ob- 
served decorrelation bandwidths greatly exceed modeled 
expectations, they are strongly diminished when scaled 
to CHIME frequencies (400-800 MHz), suggesting the 
effect would be unimportant in the CHIME band. 

A Galactic latitude dependence may also be expected 
if there are systematic errors in Milky Way electron den- 
sity distribution models (Cordes & Lazio 2002; Yao et al. 
2017). These could lead to Galactic RRATs and radio 
pulsars, predominantly a Plane population, misidenti- 
fied as extragalactic FRBs. Specifically, underestimates 
of the Galactic DM contribution in the Plane would pref- 
erentially result in higher Plane FRB rates, and the con- 
verse. That we find no latitude dependence implies any 
inaccuracies in the DM models are not strong enough to 
meaningfully impact the CHIME/FRB detection rate as 
a function of latitude. Given the subdominance of pulse- 
broadening effects (see §3.4), we caution that these lim- 
its are not as stringent as one might hope. The lack 
of latitude dependence is nevertheless interesting given 
that these models are calibrated very differently at high 
Galactic latitude versus low. For example, of the 189 in- 
dependent distance estimates used by Yao et al. (2017) 
to calibrate their Galactic DM model, for |b| < 5°, 77% 
are from either H I kinematics measurements or asso- 


ciated nebulae, while these two types are used for only 
4% of calibrators at |b| > 5°. Our results therefore sug- 
gest that either their calibration methods are as reliable 
as they believe, that any problems with methods fortu- 
itously cancel, or that the current sample size and fore- 
ground modeling for CHIME/FRB is insufficient (given 
the relative insensitivity to pulse-broadening effects). A 
lack of evidence for significant systematics in these mod- 
els suggests that efforts to use FRBs to constrain the DM 
content of the Galactic halo are worthwhile. 


6. CONCLUSIONS 


In this work we have taken a large subset of events 
from the first CHIME/FRB catalog, occurring dur- 
ing relatively uniform observing conditions, with the 
goal of determining whether a Galactic latitude depen- 
dence exists in the FRB sky distribution. We con- 
structed a radiometer-based relative sensitivity map 
over the CHIME/FRB field of view which accounts for 
sky temperature, bandwidth, intensity beam response, 
and pulse-broadening effects. With this map, we choose 
a threshold, which sets the field of view and allows us to 
filter out events considered below completeness. We per- 
formed one-sample Kolmogorov-Smirnov and Anderson- 
Darling tests over a range of sensitivity cuts, but find no 
evidence to reject the null hypothesis of isotropy with 
respect to Galactic latitude. We then partition the sky 
into high and intermediate latitude regions and perform 
a Bayesian model selection based on a biased coin flip- 
ping treatment, concluding the weight of evidence for 
isotropy ranges from substantial to strong (depending 
on the sensitivity threshold). Finally, we explored a va- 
riety of cross-checks but find no strong evidence for any 
deviation from isotropy. Our results are thus consistent 
with FRBs being a cosmological population, and we find 
no evidence for differential systematic errors in Galactic 
DM models, in spite of them having calibration strate- 
gies that depend strongly on latitude. 
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